Nature of 45° vortex lattice reorientation in tetragonal superconductors 
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The transformation of the vortex lattice in a tetragonal superconductor which consists of its 45° 
reorientation relative to the crys tal axes is studied using the nonlocal London model. It is shown that 
the reorientaton occurs as two successive second order (continuous) phase transitions. The transition 
magnetic fields are calculated for a range of parameters relevant for borocarbide superconductors in 
which the reorientation has been observed. 



Properties of the vortex matter have recently attracted 
great attention due to diversity of phases and novel phe- 
nomena associated with them. One of the main research 
goals is determination of the phase diagram. In high tem- 
perature superconductors the vortex matter phases in- 
clude the vortex liquid and various vortex solids which ex- 
ist due to the competition of intervortex interactions with 
fluctuations both thermal and those due to the quenched 
disorder.^ On the other hand, in borocarbide supercon- 
ductors a rich variety of quite perfect vortex crystals 
has been observed. The experimental information comes 
from such different measurements as neutron diffraction, 
decoration, and scanning tunneling microscopy.Bcl For 
these near isotropic materials the entropy contribution 
to the free energy is small and phase transitions in the 
vortex lattice are governed by competition between in- 
tervortex interactions of different symmetry. 

The borocarbides are materials of the tetragonal sym- 
metry. Interactions of this symmetry should exist for any 
physical subsystem of the crystal. In particular, in the 
mixed state with the field along the fourfold tetragonal 
axis they would favor a square vortex lattice. However, 
the standard magnetic repulsion of vortices is isotropic 
in this case. The isotropic interaction becomes domi- 
nant when the intervortex distance is large enough and 
a sparse lattice is close to hexagonal, the most closely 
packed two dimensional lattice. One, therefore, expects 
that the interplay of the interactions of different sym- 
metries may result in structural transformations of the 
vortex lattices, observed in borocarbides. 

For the applied magnetic field along the fourfold sym- 
metry axis, these transformations are as follows. With 
decreasing magnetic field, the lattice undergoes a sec- 
ond order phase transition, at which the square structure 
loses stability and becomes a rhombic (distorted hexag- 
onal) vortex lattice.EHa As the field further decreases, the 
rhombic lattice changes the orientation relative to under- 
lying crystal by 4EL° which has been classified as the first 
order transition.Er&Q For the field along the twofold axis, 
the 90° reorientation has been reported.0 

In this paper we study in detail the 45° reorienta- 
tion and clarify its nature. We show that this reorienta- 
tion proceeds as two successive second order (continuous) 
transitions and not as an abrupt first order transition, the 
scenario assumed before. Instead of considering a limited 



class of rhombic lattices,El we study the general class of ar- 
bitrary lattices. We find that in the field region between 
the two second order phase transitions, the lattice with 
the lowest possible symmetry is realized (with the inver- 
sion being the only symmetry element). This intermedi- 
ate region is quite narrow and the structural evolution 
in this field domain might be difficult to discern exper- 
imentally. However, the thermodynamic characteristics 
of the superconductor are different for the two scenarios, 
and this can be tested. In particular, no latent heat is ex- 
pected during the lattice reorientation. We also predict a 
peak in the critical current in the transition region if the 
pinning is of a weak collective type. Below we describe 
the London model with nonlocal corrections relevant for 
the mixed state of borocarbides. Then, the numerical 
procedure is outlined and the results are presented 

A fruitful approach to the problem of the vortex lattice 
phases is the extended London model. We start here with 
London equations corrected for nonlocalityficl 
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Here, a 3 = Aj 



($o/2tt)3j0, Aj is the vector potential, 6 
is the order parameter phase, and <f>o is the flux quantum. 
The nonlocal response kernel ?y(k) is expanded up to 
the second order terms in the wave vector k. The tensor 
fiijim {viVjViv m )"/(T,£) where v is the Fermi velocity 
and the function 7 decreases somewhat with temperature 
and drops fast for short mean-free paths £u It is difficult 
to accurately estimate the components of h because of 
uncertainties in determination of Fermi velocities and, in 
particular, of the mean-free path. At low temperatures, 
h ~ 7/K 2 , where k is the Ginzburg-Landau parameter. 
Since good crystals of borocarbides are clean materials 
with k — 10-^15, one expects the components of n to be of 
the order 10~ 2 . Note also that for the problem of vortex 
lattices in fields well under the upper critical field, the 
correction X 2 hk 2 ~ £ 2 fc 2 <C 1 (£0 is the zero-T coherence 
length). Therefore, for strong type-II superconductors, 
the corrections to the standard London equations and 
the truncation in the expansion (Q) are well justified. 

For the tetragonal symmetry, the tensor h in the crys- 
tal frame has four independent components n xxxx , n xxyy , 
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n ZZZZl and n xxzz The inverse mass tensor has two differ- 
ent components m xx = m yy and m~} . The London free 
energy functional corresponding to Eq. ([!]) reads 

F = J (\h\ 2 ~ X 2 e l j k ei mn kjk m ql k 1 h n hi^ (2) 

where h(k) is the magnetic field and e^fe is the unit an- 
tisymmetric tensor. The nonlocal corrections preserve 
linearity of the London equations and do not change the 
standard London result that the interaction of two vor- 
tices is proportional lo the field of one of them at the 
location of the other As usual, the free energy density 
of a vortex lattice is given by F — (B 2 /87r<I>o) Yl g ^z(g) 
where B is the magnetic induction, g is a vector of the 
reciprocal lattice and h z is the component of the single 
vortex field along the vortex axes. We are interested in 
the field along the fourfold symmetry axis z. Solving 
Eq. (]]]) for a singlej-vortex one can bring the free energy 
density to the formQ 

where n = n xxyy and d — 2(n xxxx — 3n xxyy ). The free 
energy F(B, T) is the thermodynamic potential, which 
is minimum in equilibrium of a superconducting slab in 
a perpendicular applied field. The temperature enters 
F(B,T) via T dependent parameters A(T), n(T) and 
d(T) that can, in principle, be calculated using a mi- 
croscopic model. Note that besides the factor B 2 , the 
induction enters via the area of the primitive lattice cell. 
We determine the stable lattice by numerical minimiza- 
tion of F(B,T;g) with respect to the lattice structure 
specified by a given set of g's. 




FIG. 1. General vortex lattice and its orientation relative 
to the crystal. 

The vortex lattice is completely defined by the basis 
vectors ai and a2, i.e., by four parameters. Since a unit 
cell accommodates one flux quantum, aia 2 sin/? = <Pq/B, 
three parameters suffice. Following Ref. [n] we choose a, 
p = (a?,/ ax) cos/3, and a = (02/01) sin (3 as the needed 



three (see Fig. 1 for definitions of a and (3). The pa- 
rameters p and a are convenient because one can select 
a domain of their variation, each point of which corre- 
sponds to a lattice with various equivalent choices of the 
basis vectors ai^. Thus, the minimization of F is done 
at fixed B, n, and d with respect to *p,a and a for 
0<a<7r, < p < 0.5, and p 2 + a 2 > 1.E3 The minima 
of F are often located on the boundaries of this domain; 
we use the "Amoeba." numerical routine convenient in 
such circumstances.Eil The cutoff factor exp (— £ 2 g 2 ) was 
introduced inside the sum (^) to properly account for the 
failure of the London model in the vortex core. Changing 
parameters B and n, d we obtain the phase diagram. 
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B, magnetic induction in units of <£>q/(2kX) 2 

FIG. 2. Phase diagram of the vortex lattice in the region 
of the reorientation for n = 0.015. Nonlocal parameters n 
and d are defined in the text. The magnetic induction b is in 
units of $o/(2ttA) 2 . 

The main finding of this work is that the reorientation 
of the lattice proceeds as two steps. Figure 2 shows the 
transition lines on the B,d plane for a fixed n = 0.015. 
The equilibrium lattices both before and after the reori- 
entation have the rhombic symmetry Z) 2 ^-li3 Their sym- 
metry axes, which coincide with the diagonals of a rhom- 
bic unit cell (with the appropriate choice of such a cell, 
see Fig. 2), are aligned with [110] and [110] at lower mag- 
netic inductions, whereas the symmetry axes are at [100] 
and [010] for higher JS's. This result is in accordance with 
data for Y N12B2CB In a narrow region between the two 
rhombic phases, a less symmetric lattice is stable. Here, 
the unit cell is a general parallelogram. All in-plane sym- 
metry elements disappear except the inversion, and the 
symmetry group reduces to Cih- 

One could describe the reorientation process as a grad- 
ual rotation of the unit cell accompanied by a slight 
deformation. Figure 3 shows how the angles a and [3 
change when B increases in the vicinity of the reorien- 
tation (for d = 0.05 and n = 0.015). The transitions 
atB« 3.18 $ /(2ttA) 2 and B as 3.27 $ /(2ttA) 2 are seen 
clearly. The angles a, (3, and the other lattice parameters 
are continuous at the two transition fields. We conclude 
that both phase transitions that occur during the reorien- 
tation are of the second order. The sequence of symmetry 
changes with the field decreasing is D^h —> C*2h — ► -D 2 /j. 
While at the first step of the reorientation, the symme- 
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try becomes lower, at the second step the symmetry in- 
creases. Correspondingly, the ground state is double de- 
generate in both D 2 h phases; the degenerate vacua (two 
equilibrium structures of the same energy) are related by 
rotations over 90°). The structure becomes four times 
degenerate in the intermediate C 2 h phase (rotations over 
±45° and 90°). Practically, this may lead to apparently 
increased disorder in the C 2 h phase. 
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B, magnetic induction in units of Oq/(2ti^) 2 

FIG. 3. Evolution of angles a and j3 (defined in Fig. I) 
with field b — <I>o/(27rA) 2 during reorientation. 



domain of monoclinic phase shrinks. Then, it would be 
difficult to distinguish experimentally this situation from 
1st order transition, because the entropy would change 
fast with B during the reorientation (for the B-sweep 
at a fixed T). Still, one should not observe hysteresis, 
characteristic of the 1st order transitions. If the sequence 
of transitions we suggest here is found, it would be of 
interest to suppress n and d by making the mean-free 
path shorter and to see how the transition evolves (as 
has beemdoiie with doping Lu-based borocarbide crystals 
with Co&M) 
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B, magnetic induction in units of Oo/(27i^) 2 



It is worth noting that the relative energy differences 
between the equilibrium C 2 h lattice and the rhombic ones 
is exceedingly small. As an example, we provide this fig- 
ure for d = 0.05 and n = 0.015: the relative difference 
between energies of the rhombic lattice at the transition 
point and the lattice in the middle of the field domain of 
C 2 h structure is of the order 10 -7 . This is much smaller 
than 10 -2 for the relative energy differences usually cited 
for triangular and square lattices within the standard 
London or Ginzburg-Landau models. 

The location of the phase transition lines is sensitive 
to both n (for the isotropic correction in Eq. (||)) and d 
(of the four-fold symmetric correction). Figure 4 shows 
the phase diagram of the vortex lattice on B, d and B, n 
planes in the region of the reorientation process. The 
region of stability of the monoclinic lattice is broader for 
smaller n's and larger d's. Still, as is seen at Fig. 4, 
this field region is narrow for values of n and d adopted 
in our simulations. For example, for LuNi 2 B 2 C with 
A w 710 A, the field unit $ /(2ttA) 2 is about 100 G. 

The new scenario of the lattice reorientation in a 
tetragonal superconductor which is found this paper has 
implications for thermodynamic characteristics of the 
vortex lattice and its dynamic behavior. We have found, 
and this is our main result, that the reorientation pro- 
ceeds as two successive phase transitions of the second 
order when the applied field or temperature vary. There- 
fore, continuous variation of the entropy (i.e., no latent 
heat) and of the reversible magnetization are expected 
during reorientation. In contrast, the old scenario of the 
first order transition implied discontinuous jumps of the 
above quantities. 

As is seen at Fig. 4, for small values of n and d, the 




1 2 3 4 5 

B, magnetic induction in units of On /(2n\) 2 

FIG. 4. Phase diagram of the reorientation transformation 
(a) in the n, 6 plane for a set of d's indicated, and (b) in the 
d, b plane for a set of n's indicated. 

Both the upper and lower phase transitions (D 2h «-> 
C 2 h) of the reorientation process cause uniform sponta- 
neous deformations of the vortex lattice. As a result, a 
particular combination of the elastic lattice moduli van- 
ishes at the transitions. It has been recently shown that a 
change of elastic properties of this type leads generally to 
peculiarities in the critical current, provided a weak col- 
lective pinning operates in the materially Therefore, the 
reorientation of the vortex lattice in borocarbides may 
lead to a peak in the critical current. 

Finally, we would like to point to other possible ap- 
plications of our results. The London model we em- 
ployed reflects properly the symmetry of the system. It 
was originally derived for an anisotropic Fermi surface 
and isotropic superconducting pairing However, the d- 
wave type of, pairing also leads to a similar effective Lon- 
don modelll3 (for not very low temperature where the 
effects of the order parameter nodes become essential). 
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The reorientation of the vortex lattice has indeed been 
found theoretically, and characterized as the first order 
transition!!- 7 ] It would be of interest to check whether or 
not our scenario of the reorientation applies to this case 
as well. 

This work is supported by NSC of Taiwan through the 
grants #89-2112-M-009-0016 and #89-2112-M-009-039. 
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